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Abstract. We review recent developments regarding radiation driven mass loss from OB-stars. We 
first summarize the fundamental theoretical predictions, and then compare these to observational 
results (including the VLT-FLAMES survey of massive stars). Especially we focus on the mass 
loss-metallicity dependence and on the so-called bi-stability jump. 

Subsequently we concentrate on two urgent problems, weak winds and wind clumping, that have 
been identified from various diagnostics and that challenge our present understanding of radiation 
driven winds. We discuss the problems of "measuring" mass-loss rates from weak winds and the 
potential of the near infrared, Br-alpha line as a tool to enable a more precise quantification, and 
comment on physical explanations for mass-loss rates that are much lower than predicted by the 
standard model. 

Wind clumping, conventionally interpreted as the consequence of a strong instability inherent 
to radiative line-driving, has severe implications for the interpretation of observational diagnostics, 
since derived mass-loss rates are usually overestimated when clumping is present but ignored in the 
analyses. Simplified techniques to account for clumping indicate overestimates by factors of 2 to 10, 
or even more. If actually true, these results would have a dramatic impact on the evolution of, and 
the feedback from, massive stars. We discuss ongoing attempts (including own work) to interprete 
the corresponding observations in terms of more sophisticated models. By allowing for porosity in 
density and velocity space, and for a non-void inter-clump medium, such models might requke only 
moderate reductions of mass-loss rates. 
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1. INTRODUCTION 

Massive stars are critical agents in galactic evolution, both in the present and in the early 
Universe (e.g., re-ionization and first enrichment). Mass loss is a key process, which 
modifies chemical profiles, surface abundances, and luminosities. Furthermore, mass 
loss has to be understood quantitatively in order to describe and predict massive star 
evolution in a correct way. The standard theory to describe hot, massive star winds is 
based on radiative line-driving, and has been proven to work successfully in most evolu- 
tionary phases (OB-stars, A-supergiants, and LBVs in their "quiet" phase). Also for the 
pivotal Wolf-Rayet (WR) stadium, line-driving is still the most promising acceleration 
mechanism 

In this review, we summarize fundamental predictions of the theory, as well as corre- 
sponding observational evidence, and subsequently concentrate on two urgent problems 
that challenge our understanding of line-driven winds, the so-called weak-wind prob- 
lem and wind clumping. We concentrate on the winds from "normal" OB-stars in all 



evolutionary phases (for corresponding results and problems regarding WR-winds and 
additional material, see the contributions by Hamann and Hillier, this volume). 



2. LINE-DRIVEN WINDS FROM HOT STARS - THEORETICAL 

PREDICTIONS 

To be efficient, radiative line-driving requires a large number of photons, i.e., a high 
luminosity L. Since L oc Tef/i?*^, not only OB-supergiants, but also hot dwarfs and A- 
supergiants undergo significant mass loss via this mechanism. Typical mass-loss rates 
are of the order of M ^ 0.1 . . . 10 -lO^^Moyr^^ with terminal velocities Vc» ~ 200 
. . . 3,000 kms^. Another prerequisite is the presence of a multitude of spectral lines, 
with high interaction probabilities, close to flux maximum, implying that the strength of 
line-driven winds should strongly depend on metallicity. ^ 

Pioneering work on this subject was performed by Lucy and Solomon [3] and Castor, 
Abbott & Klein (Hi, "CAK"), where the latter still builds the theoretical foundation 
of our present understanding. Improvements with respect to a quantitative description 
and first applications were provided by Friend and Abbott [5] and Pauldrach et al. |@], 
whereas recent reviews on the topic have been published by Kudritzki and Puis [0] and 
Puis etal. M\- 

The principle idea of radiative line-driving relies on two processes. 

1. Momentum is transferred to the wind matter via line absorption/emission processes, 
mostly resonance scattering, with a net change in radial momentum 

h 

APradial = - ( Vin COS 0in - Vout COS 0out) (1) 
C 

where Vin and Vout are the frequencies of the absorbed and emitted photons, and is 
the angle between the photon's direction and the radial unit vector. Thanks to the fore- 
aft symmetry of the emission process, on average (cos 0out) = 0, whereas (cos 0in) ~ 1, 
since (most) of the absorbed photons originate from the stellar surface. Thus, (APj-adiai) ^ 
hvin/c, and the total radiative acceleration exerted on a mass element Am per time 
interval At can be derived from considering all participating lines, 

(AP) ^ 

^ _ \''^J^/tot _ amines ^2-) 

^'^^ AtAm AtAm 

2. Due to the huge number of metallic lines as compared to the few dozens from hy- 
drogen and helium, mostly just the metal ions are directly accelerated. Their momentum 
needs to be transferred to the bulk plasma (H, He), via Coulomb collisions. The veloc- 
ity drift of the metal ions with respect to H/He is compensated for by a frictional force 
("Stokes law") as long as the ratio between drift and thermal velocity is small (e.g.. 



||9|, llO, lUj]). Otherwise (at very low wind-densities) the metallic ions might decouple 
from the wind, and the wind no longer becomes accelerated. 

The real challenge is to evaluate Eq. [21 (see also Owocki and Gayley, this volurne). 



Following CAK, this is conventionally done by (i) applying the Sobolev theory [|12|] 



to approximate the line optical depths and thus the interaction probabilities, and (ii) 
to replace the summation by appropriate integrals over the line-strength distribution 
(resulting from detailed NLTE calculations), where the line-strength k is the line-opacity 
measured in units of the Thomson-scattering opacity. This distribution can be fairly 
well approximated by a power-law, dA^(fc)/dfc oc M^ff k^zL; with A'eff the effective (flux- 



weighted) number of lines and a ^ 0.6. ..0.7 (e.g., 1)1311 ). Note that both quantities 
depend on metallicity and spectral type. As a final result, gj-ad °^ ((dv/dr)/p)", i.e., 
depends on the spatial velocity gradient and on the inverse of the density. 

Scaling relations and WLR. Once the above quantities are inserted into the hydrody- 
namic equations (adopting stationarity), the latter can be solved (almost) analytically, 
returning the following scaling relations for mass-loss rate, velocity law, and terminal 
velocity: 

M o. Ae/«V/«'(M(l-r))^"^/"', vir)=v^{l-^y 
~ 9 9^ « /2GM(l-r)x2 

Voo ^ 2.25- Vesc, Vesc= • (3) 

I — a \ / 

with Eddington-F, (photospheric) escape velocity Vesc, and a' = a — 5, where 5 ~ 0.1 
describes the run of the ionization [14]. The velocity-field exponent, /3, is of the order 
of 0.8 (for 0-stars) to 2 (for BA-supergiants). 

Using these scaling relations, a fundamental prediction for line-driven winds becomes 
apparent if one calculates the so-called modified wind-momentum rate, 

Mv^(i?,/i?0)i/2 oc A^^i/«'l1/«' (m(1 -r))'^'"'^"', (4) 

and accounts for the fact that a' is of the order of 2/3. Then the wind-momentum 
rate becomes independent on mass and F, and can be expressed in terms of the wind- 
momentum luminosity relation (WLR), discovered first by Kudritzki et al. [[l5ll . 

log(^Mv^(i?,/i?0)i/2) ^ .^;log(L/LQ) spectral type) (5) 

with slope X — \/ a' and offset D, which depends on A'eff and thus on metallicity z and 
spectral type. Originally, it was proposed to exploit the WLR for measuring extragalactic 
distances on intermediate scales (up to the Virgo cluster), but nowadays the relation is 
mostly used to test the theory itself (see below). 

Theoretical 1-D models. Though the basic scaling relations for line-driven winds are 
known since the key paper by CAK (and updates by 1 14, 1,0]), quantitative predictions 



require consistent NLTE/radiative-transfer calculations, to derive the line-force as a 
function of spectral type and metallicity, as well as the inclusion of processes neglected 
in the original work, for example line-overlap (e.g., [16, , 17]). 

The most frequently cited theoretical wind models (stationary, 1-D, homogeneous) 
are those from Vink et al. [Tsi T^. Based on the Monte-Carlo approach developed by 



Abbott and Lucy [|20D, they allow multi-line effects to be considered. In these models. 



the mass-loss rate is derived (iterated) from global energy conservation, whilst the (/3-) 
velocity field is pre-described and the NLTE rate equations are treated in a simplified 



way. Pauldrach [21] and Pauldrach et al. [|22|j23l]. on the other hand, obtain a consistent 



hydrodynamic solution by integrating the (modified) CAK equations based on a rigorous 
NLTE line-force using Sobolev line transfer. Moreover, Krticka and Kubat [[lol l24, 25 



and Krticka [27] solve the equation of motion by means of a NLTE, Sobolev line- 
force, including a more-component description of the fluid (accelerated metal ions plus 
H/He) that allows them to consider questions regarding drift- velocities, non-thermal 
heating, and ion decoupling. Also, Kudritzki [28] (see also Kudritzki et al. [29]) provides 
an analytic "cooking recipe" for mass-loss rate and terminal velocity, based on an 
approximate NLTE treatment, and Grafener and Hamann [1, 2] obtain self-consistent 
solutions (applied to WR winds) by means of a NLTE line-force evaluated in the 
comoying frame (see Mihalas et al. [3u. [3lLl32Ll33|] and Hamann, this volume). Finally, 



Lucy lU 35] and Miiller and Vink [36']derive the wind-properties from a regularity 
condition at the sonic point, in contrast to most other solutions that invoke a singularity 
condition at the CAK-critical point of the wind (see comments by Owocki, this volume). 

Results and predictions from hydrodynamic modeling. Most of the various approaches 
yiel d consistent results, e.g., when comparing the "mass-loss recipe" from Vink et al. 
[18] with similar investigations utilizing different codes [28i |23i ^3]. Moreover, the 



WLR concept is impressively confirmed by the simulations performed by Vink et al.: 
The obtained modified wind-momenta follow an almost perfect power-law with respect 
to stellar luminosity alone, independent of luminosity class, and, for solar abundances, 
"only" two distinct relations covering the complete spectral range have been found, one 
for 50 kK > Teff > 27.5 kK and the other for 22.5 kK > T^ff > 12 kK, respectively. In 
other words, the spectral type dependence of x and D in Eq. [5] seems to be rather mild. 

Also regarding the predicted metallicity dependence, the various results agree satis- 
factorily (note that the z-dependence of Vcx, is rather weak): 



Kudritzki S: °c ^0.12^ Krticka [(270 : Voo 



Vink et al. M oc z^'^'^ for 0-stars, M oc z^'^"^ for B-supergiants, 



Krticka [S: M oc z*^-^^ for 0-stars. 



3. OBSERVATIONS VS. THEORY 

In the last decade, various spectroscopic NLTE analyses of hot stars and their winds have 
been undertaken, in the Galaxy and in the Magellanic Clouds, in the UV, in the optical, 
and in a combination of both. For a compilation of these publications (without Galactic 
Center objects and objects analyzed within the FLAMES survey of massive stars, see 
below), see Tables 2 and 3 in Puis nS m. to be augmented by the UV-Pv investigation of 
Galactic 0-stars by FuUerton et al. [|38|], the UV-i-optical analysis of Galactic 0-dwarfs 
by Marcolino et al. [39], and the optical analysis of LMC/SMC 0-stars by Massey et al. 
[I4O]. Most of this work has been performed by means of 1-D, line-blanketed, NLTE, 
atmosphere/spectrum- synthesis codes allowing for the presence of winds, in particular 
CMFGEN (Hillier and Miller ^^), WM-Basic (Pauldrach et al. ^\), and FASTWIND 
(Puis et al. (4^). 



Central results. The results of these investigations can be roughly summarized as fol- 
lows, (i) The mass-loss rates from SMC stars (with z ~ 0.2 zq, see [43,1 and refer- 
ences therein) are indeed lower than those from their Galactic counterparts, (ii) For 



O- and early B-stars, the theoretically predicted WLR from Vink et al. [|18|1 is met, ex- 
cept for 0-supergiants with rather dense winds, in which the observed wind-momenta 
are higher (by factors around three) than the predictions (which might be explained by 
wind-clumping effects, see Sect. [S]), and for a number of late 0-dwarfs (and a few O- 
giants), in which the observed wind- momenta are much lower than the predictions (this 
is the so-called "weak-wind problem", see Sect. Hj). (iii) B-supergiants below the "bi- 
stability jump" (Tgff < 22 kK) show lower wind-momenta than predicted, as outlined in 
the following. 

The bi-stability jump: predictions and observations. A fundamental prediction by Vink 



et al. [1180 is the occurrence of two distinct WLRs, one for hotter objects and one for 
cooler objects, with the division located around 25±2.5 kK. This rather abrupt change is 
due to the so-called bi-stability mechanisrrQ, which relies on the fact that the mass-loss 
rates of line-driven winds are, for typical chemical compositions, primarily controlled 
by the number and distribution of /ron-lines, because of their dominant contribution 
(~50%) to the total line acceleration in the lower wind [fl3l 18 . l27ll . Below roughly 



25 kK, the ionization of iron is predicted to switch abruptly from Fe IV to Fe III, and 
since Fe ill has more driving lines than Fe iv at flux maximum, the mass-loss rate must 
increase. Quantitatively, Vink et al. [18] predict an increase in M by a factor of five and 
a decrease of Voo by a factor of two, so that, overall, B-supergiants (except for the earliest 
sub-types) should have higher wind-momenta than their 0-star counterparts at the same 
luminosity. 

Observations confirm the "velocity-part" of this picture, at least qualitatively. For stars 
with reff> 23 kK, the observed ratio is Voo/vesc ~ 3, whereas it decreases gradually 
towards cooler temperatures, reaching values of Vcx,/vesc~ 1.3. .. 1.5 for stars with reff< 



18 kK [|45L l46L |47|1 . With respect to the predicted increase in M, however, the situation 



is different. As shown by Markova and Puis [|47[1 . the mass-loss rates of B-supergiants 
below the observed bi-stability jump (reff< 22 kK) actually decrease or at least do no 
change. This is a first indication that there are still problems in our understanding of 
line-driven winds. 

The FLAMES survey of massive stars. Further progress has been obtained within the 
FLAMES survey of massive stars (P.I. S. Smartt), a project that performed high resolution 
multi-object spectroscopy of stars located within eight young and old clusters in the 
Galaxy and the Magellanic Clouds. In total, 86 0-stars and 615 B-stars were observed 
(for introductory papers and a brief summary, see Evans et al. [48i,i49i,^50]). The major 
scientific objectives of this survey were to investigate (i) the relation between stellar 
rotation and abundances (i.e., to test the present theory of rotational mixing), (ii) the role 
of binarity, and (iii) stellar mass-loss as a function of metallicity. 

Regarding the last objective, Mokiem et al. lIsTl 152 ] analyzed a total of ~ 60 O- 
and early B-stars in the SMC and LMC, by means of fastwind and using a genetic 



denoted after some peculiar behaviour of theoretical models for the wind of P Cygni |[44ll . 



algorithm 115 311 . The resuks were combined by Mokiem et al. [|43ll with data from 
previous investigations, to infer the metallicity dependence of line-driven mass-loss 
based on a significant sample of stars. Using mean abundances of z = 0.5 zq (LMC) 
and z = 0.2 zq (SMC), a metallicity dependence of oc (z/zq)^'^^, and a correction for 
clumping effects (see below) following Repolust et al. [54], they derived an empirical 
relation 



M (z/zq) 



0.72±0.15 



(6) 



with rather narrow confidence intervals. This result is consistent with theoretical predic- 
tions, both from line-statistics [Il3n and from hydrodynamic models (see above). 



4. WEAK WINDS 

The results as summarized above imply that line-driven mass loss seems to be basically 
understood, though certain problems need further consideration. In particular, from early 
on there were indications that the (simple) theory might break down for low-density 
winds. E.g., Chlebowski and Garmany i55ll have derived mass-loss rates for late O- 
dwarfs that are factors of ten lower than expected. By means of UV-line diagnostics, 
Kudritzki et al. [56] and Drew et al. [57] have derived mass-loss rates for two BII stars 
that are a factor of five lower than predicted, and Puis et al. [|58l] have shown that the 
wind-momentum rates for low-luminosity dwarfs and giants (logL/L© < 5.3) lie well 
below the empirical relation for "normal" 0-stars. 

The last investigation illuminated an immediate problem arising for low-density 
winds. For M < (5. . .1) ■ 10^^ Mqjv^^, the conventional mass-loss indicator, H^, be- 
comes insensitive, and only upper limits for M can be derived (for a recent illustration 
of this problem, see [39]). Instead, unsaturated UV resonance lines (Civ, Siiv, Cm) 
might be used to obtain actual values for M (e.g., f,59, ,8„.39]). 

By means of such UV-diagnostics, strong evidence has accumulated that a large 
number of late type 0-dwarfs (and a few giants of intermediate spectral type) have 
mass-loss rates that are factors of 10 to 100 lower than corresponding rates from both 
predictions and extrapolations of empirical WLRs. In particular, such weak winds have 
been found in the Magellanic Clouds (0-dwarfs in NCG346 (LMC): Bouret et al. @; 
extremely young 0-dwarfs in N81 (SMC): Martins et al. [|59|]) and in the Milky Way 
(0-dwarfs and giants: Martins et al. (6]^; late 0-dwarfs: Marcolino et al. ^39]). 

Two points have to be stressed, (i) Until now, it is not clear whether all or only part of 
the late type dwarfs are affected by this problem, (ii) The derived UV mass-loss rates are 
not very well constrained, since they might be contaminate^l from X-rays embedded in 
the wind (due to shocks, see next Section). The higher the X-ray emission, the weaker 
the lines, and the higher the actual mass-loss rates (see Figs. 19 and 20 in Puis et al. [S]]). 
However, to "unify" the present, very low, M-values with "normal" mass-loss rates by 
invoking X-rays, unrealistically high X-ray luminosities would be required [|39il . 

The weak-wind problem is a prime challenge for the radiative line-driven wind theory. 
Martins et al. [|59ll investigated a variety of candidate processes (e.g., ionic decoupling. 



via a modified ionization equilibrium. 



shadowing by photospheric lines, curvature effects of velocity fields), but none of those 
turned out to be strong enough to explain the very low mass-loss rates that seem to be 
present. At the end of this review, we will return to this problem. 



5. WIND CLUMPING 

During the last years, overwhelming direct and indirect evidence has accumulated that 
one of the standard assumptions of conventional wind models, homogeneity, needs to 
be relaxed. Nowadays the winds are thought to be clumpy, consisting of small scale 
density inhomogeneities, where the wind matter is compressed into over-dense clumps, 
separated by an (almost) void inter-clump medium (ICM). Details on observations and 
theory can be found in the proceedings of a recent workshop, "Clumping in hot star 
winds" [62]. 

Theoretically, such inhomogeneities are considered related to structure formation due 
to the line-driven ("de- shadowing") instability, a strong instability inherent to radiative 
line-driving (cf. Owocki, this volume). Time-dependent hydrodynamic models allow- 
in g fo r this instability to operate have been developed by Owocki and coworkers (1- 
D-H HQ; 2-D:H |^) and by Feldmeier H and show that the wind, for 



r ^ l.3R^, develops extensive structure consisting of strong reverse shocks separating 
slower, dense material from high-speed rarefied regions in between. Such structure is 
the most prominent and robust result from time-dependent modeling, and the basis for 
our interpretation and description of wind clumping. Within the shocks, the material 
is heated to a couple of million Kelvin, and subsequently cooled by X-ray emission 
(which has been observed by all X-ray observatories), with typical X-ray luminosities 
Lx/Lboi ~ 10^^ (for newest results, see Sana et al. [70]). 

Clumping effects. Until now, most diagnostic methods to investigate the effects of clump- 
ing use the following assumptions: The clumps are optically thin, the ICM is void, the 
velocity field remains undisturbed, and the so-called clumping factor, /d, measures the 
over-density inside the clumps with respect to the average density. This simple model of 
micro-clumping allows one to incorporate clumping into NLTE-codes without any major 
effort, namely by multiplying the average (wind-) density by /d and by multiplying all 
opacities/emissivities by the inverse of /d (i.e., by the volume filling factor). 

The most important consequence of such optically thin clumps is a reduction of any 
M derived from p^-dependent diagnostics (e.g., recombination based processes such 
as Wa or radio-emission), assuming smooth models, by a factor of V7ci- That there 
is a reduction is conceivable, since, under the assumptions made, the square of the 
over-density "wins" against the smaller absorbing/emitting volume. Thus, a lower M 
is sufficient to produce the same optical depths/emission measures as in smooth models. 

Note, however, that in this scenario any M derived from p-dependent diagnostics (e.g., 
UV-resonance lines) remains uncontaminated, since in this case the over-density cancels 
against the smaller absorbing/emitting volume. Finally, it should be mentioned that a 
clumpy medium also affects the ionization equilibrium, due to enhanced recombination 



(e.g., m). 

Results from NLTE- spectroscopy allowing for micro-clumped winds are as follows, 
(i) Typical clumping factors are /ci~ 10. . .50, and clumping starts at or close to the 



wind base, the latter in conflict with theoretical predictions. Derived mass-loss rates are 
factors of 3 to 7 lower than previously thought [T^, Izli [63, [vi]. In strong winds, the 



inner region is more clumped than the outer one (/™ ?a 4. . .6 x /°" ), and the minimum 



reduction of smooth mass-loss rates is by factors between 2 and 3 117411 



The P V problem. From a mass-loss analysis using the FUV P v resonance lineH for a 
large sample of 0-stars, FuUerton et al. 113 811 (see also 117511 ') concluded that the resulting 



mass-loss rates are a factor of 10 or more lower than derived from Hce and/or radio 
emission using homogeneous models, implying 100! Similar results have been 
found from unsaturated P Cygni lines from lower luminosity B-supergiants [76|]. 

If such large reductions in M were true, the consequences for stellar evolution and 
feed-back would be enormous. Note that an "allowed" reduction from evolutionary 



constraints is at most by a factor of 2 to 4 (Hirschi Tim . 



Porosity and vorosity. A possible resolution of this dilemma might be provided by 



considering the porosity [|78|1 of the medium, also suggested to explain the observed 



X-ray line emission (cf. the contributions by Oskinova et al. and Cohen et al. in [|62l1 . 
and particularly the discussion on X-rays). Whenever the clumps become optically thick 
for certain processes, as might be true for the Pv-line, the geometrical distribution 
of the clumps becomes important (size vs. separation, shape). In this macro-clumping 
approach (see also Hamann, this volume), the effective opacity becomes reduced, i.e., 
the wind becomes more transparent ("porous"), because radiation can propagate through 
the "holes" provided by the ICM. Additionally, clumps hidden behind other clumps 
become ineffective because the first clump is already optically thick. 

Oskinova et al. [j79|l used a simple, quasi-analytic treatment of macro-clumping (still 
assuming a smooth velocity law) to investigate Pv in parallel with from C, Pup. 
Whereas macro-clumping had almost no effect on H^, since the transition is optically 
thin in the clumps, P v turned out to be severely affected. Thus, only a moderate reduc- 
tion of the smooth mass-loss rate (factors 2 to 3) was necessary to fit the observations, 
consistent with the evolutionary constraints from above. 

This model has been criticized by Owocki [80], who pointed out that not only the 
distribution/optical thickness of the clumps is important, but also the distribution of 
the velocity field, since the interaction between photons and lines is controlled by the 
Doppler-effect. Also the "holes" in velocity space, due to the non-monotonic character 
of the velocity field, lead to an increased escape (thus, he called this process velocity- 
porosity = "vorosity"), whilst the different velocity gradients inside the clumps lead to 
an additional modification of the optical depth. 

Resonance line formation with porosity and vorosity. To clarify in how far the above ar- 
guments/simulations depend on the various assumptions, and to characterize/quantify 
the various effects from inhomogeneous winds of different sub-structures, a current 
project in our group investigates the resonance-line formation in such winds. To this 
end, pseudo 2-D hydrodynamic models (based on different snapshots of corresponding 
1-D models from Owocki and Feldmeier, aligned as independent slices of opening angle 
0), as well as 2-D models based on a stochastic description, have been created (Fig. [H 
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FIGURE 1. Left: Density contours of stochastic (upper) and pseudo 2-D hydrodynamic wind models as 
investigated by our group. Right: Line profiles for an intermediate strong line formed in inhomogeneous 
winds with different sub-structures. See text. 



left panel). For these models then, a detailed Monte-Carlo line transfer (discarding the 
Sobolev-approximation) is performed. The right panel of Fig. [H shows prototypical pro- 
files from such simulations, based on the stochastic 2-D wind description, for an inter- 
mediate strong line that would be marginally saturated in smooth models (dashed). The 
grey dashed-dotted profile displays the effects of porosity alone (i.e., a smooth veloc- 
ity field has been used), using a rather low clumping factor, = 3.3, and an average 
separation of clumps ~ in the outer wind. Already here, a strong de-saturation of the 
profile is visible. The grey dashed-dotted-dotted line displays the other extreme, namely 
vorosity alone (i.e., now the density is smooth), using a stochastic description of the ve- 
locity field, characterized by a "velocity clumping factor" (as defined in [80], Fig. 1) /vei 
= 0.3. Interestingly, the de-saturation of the profile is similar to the porosity-effect alone. 
The solid black line displays the combined effect from porosity and vorosity, with a fur- 
ther de- saturation. If compared to a line from a smooth model of similar profile strength 
(dotted), it turns out that the effective opacity in the structured model(s) has been re- 
duced by a factor of 20, i.e., the actual M would be a factor of 20 higher than derived 
from a smooth model. Thus, structured models invoking porosity and vorosity might 
indeed resolve the discordance between the results by FuUerton et al. and evolutionary 
constraints. 

We note, however, that the profile- strength reduction presented in Fig. [T] corresponds 
to a "most favourable case", using rather ideal parameters. Our investigations have 
shown how details on porosity, vorosity, and the ICM, all are important for the forma- 
tion of the line profiles. In fact, the strengths of similar profiles calculated from our 
pseudo 2-D hydrodynamic models are only reduced by ^ 10%, because of insufficient 
vorosity inherent to structures from present time-dependent modeling (see also Owocki 
llsoll ). Such a modest reduction is much lower than needed to alleviate the discrepancy 
discussed above. Also, as it turns out, the ICM is a crucial parameter if to de-saturate 
intermediate strong lines and, at the same time, allowing the formation of the observed 



saturated profiles. Tests have shown that, with a void ICM, the formation of saturated pro- 
files is only possible if the average clump separation (controlling the porosity) is very 
small, but then the de-saturation of intermediate strong lines becomes marginal. Only 
by assuming an ICM with sufficient density O.Olpsmooth) we have been able to form 
saturated lines in parallel with de-saturated ones of intermediate strength. This finding 
is consistent with results from Zsargo et al. [81], who pointed out that the ICM is crucial 
for the formation of highly ionized species such as O vi. 

Further details and results from our investigations will be given in a forthcoming paper 
(Sundqvist et al., in prep for A&A), including a systematic investigation of different key 
parameters and effects. Future plans include a comparison with emission lines, and the 
development of simplified approaches to incorporate porosity /vorosity effects into NLTE 
models. 

6. WEAK WINDS AGAIN - BRa AS A DIAGNOSTIC TOOL 

In the preceeding paragraphs, we have argued that (i) mass-loss rates from unsaturated 
UV line-profiles are much lower than those from or radio emission, and that (ii) 
this discordance might be mitigated by porosity /vorosity effects. Recall here that the 
mass-loss rates from weak winds discussed so far (Sect. ID) rely on the same UV-line 
diagnostics, and the question arises whether one encounters a similar problem, i.e., an 
under-estimation of the "true" mass-loss rates due to insufficient physics accounted for 
in the diagnostics. Thus, to clarify in how far the weak wind problem is a real one, 
independent diagnostics are required! 

Already in 1969 Auer and Mihalas [82], based on their first generation of NLTE, 
hot-star model atmospheres, predicted that the IR Br^-line should show significant 
photospheric core emission, due to an under-population of its lower level {n = 4) relative 
to the upper one (n = 5), resulting from a very efficient decay channel 4 — > 3. Indeed, 
such core emission has meanwhile been observed in various weak wind candidates such 
as T Sco (B0.2V), HD 36861 (08111(f)), and HD 37468 (09.5V) (Najarro, Hanson and 
Puis, in prep, for A&A). Recent simulations (Puis et al. jsl]. Figs. 21/22) actually show 
that such photospheric -i- wind emission can fit the observations quite nicely, and that the 
core of Br^ is a perfect tracer for the wind density also for thinner winds (as opposed 
to Hce). Astonishingly, the height of the peak increases for decreasing M, which is 
related to the onset of the wind, i.e., the density /velocity structure in the transition 
zone between photosphere and wind, and not due to radiative transfer effects. The 
higher the wind-density, the deeper (with respect to optical depth) this onset, which 
subsequently suppresses the relative under-population of n = 4 due to efficient pumping 
from the hydrogen ground-state. Moreover, Br« is only weakly affected by the presence 
of X-rays, and thus an ideal tool to infer very low mass-loss rates. From fits to the 
observations, it turns out that M is actually very low (of the order of lO^^'^MQyr^^ for 
HD 37468, and even lower, if the wind-base were clumped). Thus, weak winds seem to 
be a reality! 

What may then be the origin of weak winds? Krticka and Kubat [[26|] argue that weak- 
winded stars display enhanced X-ray emission, maybe related to extended cooling zones 
because of the low wind density. Already Drew et al. [57] pointed out that strong X- 



ray emission can lead to a reduced line acceleration, because of a modified ionization 
equilibrium, and since higher ions have fewer lines. Thus, weak-winded stars might be 
the result of strong X-ray emission. Let us now speculate whether such strong emission 
might be related to magnetic fields. Note that weak winds can be strongly affected by 
relatively weak 5-fields, of the order of 40 Gauss according to the scaling relations 
provided by ud-Doula and Owocki [83], which is below the present detection threshold. 
In this case then, colliding loops might be generated, which in turn generate strong and 
hard X-ray emission in the lower wind, which finally might influence the ionization and 
thus radiative driving. Future simulations coupling magneto-radiation-hydrodynamic 
wind codes with a self-consistent description of the line-acceleration will tell whether 
this mechanism might work. 
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